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SUMMARY 


This report describes the first two years of work on 
a project aimed at learning how the peak and minimum pool 
boiling heat fluxes are influenced by gravity, size, and con- 
figuration, as well as other variables. 

The first months of this endeavor were devoted to the 
design and fabrication of a centrifuge facility which allowed 
a 100 _ fold variation of gravity. This apparatus is described 
briefly and the following experiments are reported: 

1. ) The Taylor unstable wavelengths are measured 
during film boiling on horizontal cylinders and their range 
is shown to be describable by adaptations of existing theory 
over a wide range of conditions , but a low gravity (or small 
size) limitation is imposed. The minimum heat flux based 
upon this wavelength is also discussed. These microstudies 
of the transition from film to nucleate boiling provide the 
basis for acceptance of the hydrodynamic theory and put some 
limitations- on the range of acceptance. 

2. ) The peak heat flux on horizontal cylinders is pre 
dieted on the basis of the hydrodynamic theory and verified 
with almost a thousand data points over a very wide variety 
of conditions. The low gravity (small size) limitation is 
corroborated in the peak heat flux case. 

3. ) About 900 observations of the peak heat flux on 
horizontal ribbons verify that induced convection can serious 
ly influence the peak heat flux. A method is developed for 
correlating this effect. 

Throughout these individual studies a general function 
al expression is developed for correlating the peak and min- 
imum heat fluxes. This relation gives the ratio of the ex- 
treme heat flux in a particular configuration, to the corres- 
ponding value on a flat plate, as function of two scale para- 
meters . Contact angle and the liquid-vapor density ratio are 
additional weak influences that can appear in the function. 

A number of observations as to the characteristic behavior 
and limitations of this expression are advanced. 


I. 


INTRODUCTION 


w 


This report summarizes a two-year study of the peak and 
minimum pool boiling heat fluxes, q m ax anc ^ ^min/ under a wide 
variety of conditions . The study was motivated by the fact that 
gravity, g, and size exert inter-related effects upon the extreme 
heat fluxes. This notion was suggested in 1963 by Costello and 
Adams [1]1 and later developed more quantitatively in references 
[2] and [3] . 


Actually, the intrinsic complexity of the effect of gravity 
began to come to light about the turn of this decade when the 
first experimental investigations of the extreme heat fluxes were 
made. Many investigators, notably Costello and Adams [4] , Siegel 
and Howell [5] and Merte et al. [6] , [7] have since measured qmax 

and q m i n under conditions of both elevated and reduced gravity. 
Their data indicate that the "one-quarter power" dependence of 
qmax and qmin upon g, that is suggested by Zuber's [8] , [9] equa- 
tions , is only approximate in some cases and inaccurate in others. 


Zuber's equations were developed upon consideration of the 
hydrodynamic stability of vapor removal during the boiling of an 
inviscid liquid on an infinite horizontal flat plate. Designat- 
ing the flat plate configuration with the“subscript , F, we can 
write Zuber ' i 

1/2 


s equations in the f orm^ : 

qmax F = 2I Pg V 2 h fg Ia g ( Pf-Pg^ 3 / 4 |( 1 + ^f) (* 


7T 


and 


qmin-p 


T T^ 
~60 


1/4 


(-) 


P h^_ 

g f g 


0 

gg (P £ ~Pg) 

. < p f + V 2 - 


16 -TT 



1 1/4 


( 1 ) 


( 2 ) 


The factor in braces in equation (1) characterizes the inflow of 
makeup liquid to compensate vapor removal. It varies from very 
near unity over most of the range of pressure, up to 1.14 at the 
critical point. Thus we can write with very good accuracy: 


g maxp 


p g 1/2 h fg [ag(p f -p g )] 1/4 


(3) 


Unless noted otherwise, qmaxp will be used to designate the est- 
imate of the peak heat flux on a flat plate given by equation 
(3). 


^Numbers in square brackets denote entries in the References 
section. 

2 Symbols not defined in context are explained in the Nomenclature 
section. 


* 


In 1963 Lienhard and Wong [10] provided some explanation 
for the failure of the relation, 

^max or qmin w % 9 1/4 , (4) 

F F 


suggested by equations (1) and (2) or (3) . They presented a hy~ 
drodynamic derivation of q m i n for a configuration with a charac- 
teristic length, namely infinite horizontal cylinders of radius, 
R. They obtained the relation 


Pcr h fcr 
= .057 g g 
R 


a min 


~ (p f -p ) 

2g — - — 2- + 


a 


iV 2 


g(Pf-Pg) 1 


[-V4 


(p f +p g ) (p f +p g )R" 


2R* 


(5) 


which differs greatly from equation (2) , especially as q m i n de- 
pends on g. 


Two years later Lienhard and Watanabe [2] noted that equa- 
tion (5) could be put in a dimensionless form which can be cast 
as follows: 


%nin 


Pfminp 


r i 

18 

|_3 . 8 

R' 2 (2R ,2 +1)_ 


1/4 


( 6 ) 


where R' is a dimensionless radius defined as 


R* s Rtgtpf-PgVa] 1 ^ (7) 

They then provided experimental data to show that, at least for 
horizontal cylinders , both 


and = f(R') (8) 

"maxp PIminF 

were true. Reference [3] went on to validate this suggestion 
using variable gravity data from a variety of sources , and in 
additional geometric configurations . In non-cylindrical geome- 
tries, the function, f(R') was replaced by f(L') where L'was 

L* = L[g(p f -p g )/0] 1 /2 (7a) 

and L was the characteristic length appropriate to the geometry 
under consideration. 

It was clear at this point that equation (4) could only be 
true for the very restrictive case in which f(L') = constant. 

The influence of gravity appeared both as a gt/4 factor in q ma x-c. 
and in the function, f(L')/ of the scale parameter, L', as wellT- 

The present work was initiated in the Fall of 1967 with the 
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following objectives in mind: 

-w 

1.) The applicability and limitations of equation (8) were 
to be fully explored. 


2 • ) Imax/ Smaxp a ? 1<3 - ^min/ Iminp versus L ' 
be developed for a variety of configurations , 
be available for design use. 


surfaces were to 
so that they would 


3 . ) The dynamics of vapor removal were to be studied in an 
effort to create additional explicit predictions of q max and/or 
q m in comparable with equations (1) , (2) , and (5) . 

Chapters IV, V, and VI discuss three experiments and anal- 
yses, each of which will bear upon all three objectives. Chapter 
II will briefly describe the variable gravity apparatus that was 
developed here for this work, and will introduce the experiments. 
Chapter III explains in greater detail our i. priori thoughts 
about the first objective, above. 


I 


II. EXPERIMENT 


The basic experimental requirement in this work was an ap- 
paratus that would permit variation of the independent variable, 
R' or L' , in equation (8) and observation of the resulting de- 
pendent variable, qmax/'Imaxp • Three factors influence L' : the 

size, L; the gravity, g; and the liquid being boiled. Thus, 
provision had to be made for mounting heaters of various sizes 
in each configuration of interest, for varying the gravity over 
a wide range, and for accomodating diverse saturated liquids. 
Provision for supplying electric power to the heater, and for 
reading this power at burnout, was needed for the dependent var- 
iable. 


The answer to this requirement was a general purpose cen- . 
trifuge apparatus which is described in detail in [11] . Figure 1 
is a general view of this facility and Fig. 2 provides a schem- 
atic view showing the various slip rings and other access on the 
underside . 

The centrifuge provides a capability for varying gravity 
between about 5 and 100 times earth normal gravity. Slip rings 
provide 2-50 amp power circuits, 3-20 amp power circuits, 
and 6 low current instrumentation circuits. A main rotating 
seal permits evacuation of a test capsule to any desired reduced 
pressure, and a secondary rotating seal facilitates measurement 
of this pressure. 

Users of the centrifuge are protected by a 0.32 cm (1/8 in.) 
thick steel cylinder from the danger of any mishap that might 
occur during operation. The foreground of Fig. 1 is designated 
the "service side" and all handling of the test section is done 
through an access door on this side. The opposite side is 
covered with plexiglas to protect the experimenter during visual 
observation of the tests. This side is designated the "observa- 
tion side" of the facility. 

Three kinds of visual studies are made. Figures 1 and 2 show 
the centrifuge rigged for one of these, namely high speed motion 
picture photography which is accomplished through two removable 
mirrors. Still photography can also be done with a camera located 
directly above the path of the centrifuge arm on the observation 
side, after removing the motion picture mirrors and light. The 
still photos are illuminated with a photoelectrical ly triggered 
strobe light directly below the camera. Visual observations were 
made by looking directly through the test section into the diffused 
strobe light. By triggering the light on each revolution of the 
arm it was possible to create the illusion of a static boiling 
process under high gravity, and to observe the transitions as 
though they were standing still. 
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Fig. 2 Schematic view of centrifuge facility 
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The test capsule used in all of the present experiments is 
shown in Fig. 3. The inside of the^ capsule is 17.8 cm (7 in.) 
long, by 7.6 cm (3 in.) high, by 9.5 cm (3.75 in.) wide. The 
tab shown in the center is cut to exactly 2.540 ± .003 cm 
(1.000 ± .001 in.) across the outer edges to provide a scale in 
photographic observations . 



Fig. 3 The test capsule mounted with a horizontal wire heater, 
and an alternate ribbon heater 


The capsule is shown loaded with a horizontal wire heater, 
10.16 cm (4.00 in.) in length — the standard length that we used 
for all of our wire heaters. The mounting schemes that we used 
for these heaters is shown in Fig. 4. These configurations were 
devised to minimize possible vapor pockets, heat loss, and 
electrical resistance. 

An alternate test heater configuration is also shown in 
Fig. 3. This is a long horizontal ribbon of limited width. It 
is 10.16 cm (4.00 in.) in length and mounted on a wide insulat- 
ing plate. 

All tests reported here were made on either of these 
heaters. The wires used in the horizontal cylinder configura- 
tion ranged from 36 to 12 gage in size. The flat ribbons ranged 
from .117 to 2.54 cm in width. 

The results of three basic experiments made on these heat- 
ers will be included in this report in support of our predictions 
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electrode 


Solder 


it 





Heater wire 


Alternate mounting for heavy gage wires 


Brass electrode 



FIG. 4 CONFIGURATION OF CYLINDRICAL HEATER WIRE ATTACHMENT 


and methods of correlation. They are: 


1 . ) Measurements of the wavelengths during film boiling on 
horizontal cylinders . These experiments are described in full 
detail in [111 and the results are reported in [12] . 

2. ) Measurements of ^max on horizontal cylinders . These 
experiments and the results are described in full detail in [13] . 

3 . ) Measurements of Imax on ribbons . These experiments are 
described in full detail in [14] and the results are reported in 
[15] . 


In general, these tests were made in five liquids: water, 
methanol, isopropanol, acetone, and benzene. Before each observ- 
ation or small set of observations, the- heater was carefully 
washed in soap and very hot water, and then rinsed in acetone. 

It was then installed in the capsule; the capsule was charged 
with liquid; and the liquid was brought to saturation and boiled 
for at least 10 minutes prior to making measurements. 
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Ill . SCALING LAWS 


The Relevant Variables 


Equation (8) is suggested by a theory which ignores viscos- 
ity and incorporates other simplifications that are not true in 
general. To obtain a broader scaling law we might first consider 
the relationship 


Smax or q min = f(h fg , p f -p g , p g , g, a, L, u, 0 C ) (9) 

In this expression, we adopt the point of view that the 
density difference (Pf-Pg) is a more fundamental variable than 
the liquid density alone. This is because it is basic to the 
buoyant force g (pf-p g ) which is important in most boiling pro- 
cesses. Certain other variables have been omitted since there 
is a good deal of experience that indicates they do not influence 
q ma x and <2 m in appreciably. These include thermal conductivities 
and heat capacities, viscosity of the vapor phase, and surface 
roughness which we shall discuss below. 

Equation (8) was based on the implicit assumptions that 
the last two independent variables, y and e c , were also unimpor- 
tant, and that p a /pf << 1. Let us see how these three assump- 
tions come about y ana what they imply. 

1.) The liquid viscosity, y, is generally ignored in hydro- 
dynamic theories of boiling since it does not impose serious 
damping in interfacial waves. However, heaters of finite size 
can induce flows around themselves through the viscous drag ex- 
erted by the rising bubbles. These flows can, in turn, disrupt 
the vapor removal pattern above a heater and seriously influence 

^max or ^min* 

We shall look only for this kind of "induced convection" 
effect of viscosity and ignore its influence on the stability 
of interfacial waves . It is a secondary effect in that the 
viscous action itself is not important in the hydrodynamic 
"crisis" of transition. Instead, viscous action acts far away 
from the immediate heater surface and it serves to create a 
flow field which in turn exerts the actual influence on the 
transitions . 
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( 10 ) 


Borishanski (16] suggested the dimensionless group^ 

n = ~~*r [tf/g(p f -p g n 1/2 

y 2 g 

to correlate viscous effects. A few years before Zuber predicted 
PlmaX]?' Borishanski established the following experimental expres- 
sion for q max on fairly large horizontal disc heaters: 

^max — (1 + 30.5/N ^ ] <3 maXF (11) 

In most cases, N is sufficiently large that 30.5/N 2//5 << 1, so the 
term in brackets proved to be a very minor correction factor. We 
shall also find that the influence of N is negligible in the hori- 
zontal cylindrical geometry. Chapter VI, however, will describe a 
situation in which N exerts great importance. 

2. ) The contact angle, 0 C , is probably the correct variable 
to characterize those surface effects that do not arise from mech- 
anical roughness. It has been established by Berenson [17] and 
Lyon [18] that surface roughness is an unimportant variable, but 
that chemical effects can cause variations in q max that reach the 
order of 10 percent on fairly large heaters under extreme circum- 
stances of surface preparation. Costello and Frea [19] present 
data that indicate these effects upon qmax can b® even more pro- 
nounced, but only in smaller scale geometri-.es. The very early 
qmax ^ ata of Cichelli and Bonilla [19] indicated that very serious 
chemical deterioration of large flat plates resulted in as much as 
a 15 percent increase in q max . 

Berenson' s observations of q m i n on flat plates showed that 
surface energy effects could cause very serious upward deviations 
of q m in from the limiting hydrodynamic minimum. 

In our work we shall deal only with data obtained on care- 
fully cleaned surfaces. It appears to be true that 0 C is very 
nearly the same (about zero) for most materials in most fluids . 
Thus cleaning will minimize (but not eliminate) the influence of 
0 c — at least for q max on comparatively large heaters. 

3. ) The density of vapor, P g , apparently exerts an indepen- 
dent influence on q max when the inflow of liquid, required to bal- 
ance vapor outflow, has to be considered. The velocity of inflow, 
which can alter the point of Helmholtz instability of escaping va- 
por jets, will be appreciable only when p g begins to approach the 
order of magnitude of Pf- Pg. At almost all saturation temperature: 
of practical interest, Pf-P g is much larger than pg as we note in 

^Perhaps the factor pf in this expression would better be written 
as (pf-p g )+p g in keeping with our emphasis upon (pf-p g ) as a more 
significant quantity than pf. 
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Fig. 5. It is this influence that was reflected in the factor in 
braces in equation (1), and (at least in that instance) the re- 
sulting effect on qmax was only a 14 percent increase near the 
critical point. It is entirely possible that the ratio, p g / (Pf-Pg), 
is generally of little consequence. 



fraction of critical temperature 


Fig. 5 Temperature dependence of Pg/ (p^-pg) for water 


The Form of Correlation Equations and Predictions 

Before we proceed to develop the more general dimensionless 
form of equation (9) , let us consider the form it takes when the 
influences of viscous drag and surface energy are neglected. 
Without y and 0 C equation (9) becomes: 

<3max or Smin = f < h fg' Pf“Pg' Pg' 5' *■> < 12 > 

Equation (12) relates seven variables in the four dimensions: 
thermal energy, length, time, and force or mass. The Buckingham 
Pi theorem shows that equation (12) must be reducible to a relation 
among any 7 minus 4 or three independent dimensionless groups of 
these variables . 

The three groups that we shall construct are: the indepen- 

dent variable, L’ , the parameter, Pg(PfPg)/ and the dependent 
variable, J 

qmax/P g 1/2 h fg [ug (p f -p g ) ] 1/4 

for q m i n it is advantageous to multiply this group by a function 
of ’p g /(p f — p g ) equal to [ (p f -p g ) p g ] V2 [1+ 2 p g / (p f -p g ) ] 1 /2 f to get: 
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qmin/Pghfglag (p f -pg)/ (pf+Pq) 2 3 1/4 

Finally, if the appropriate constants are incorporated in 
these dependent variables , the correlating equations become : 


Imax 1m in 

Imaxp or Iminp “ f[L ' Pg/ (pf"Pg) 3 


(13) 


We should like to have q m ax F and Iminp represent as closely as 
possible the actual extreme^heat fluxes on infinite flat plates. 
The dependent variable will then provide a physically meaningful 
comparison of observed values with an apt reference. There are 
few data available with which to check Zuber's equations in the 
geometry for which they were intended. However, those data 
available for heaters which approach flat plates , strongly sug- 
gest that equation (1) or (3) for q ma x is the best estimate we 
can presently make. 


Berenson showed that this is not the case for q^inf though. 
His analysis revealed that the constant, (tt2/60) (4/3)i/4, was 
incorrectly derived and had to be replaced with an experimental 
number. For this he obtained .09, so we shall henceforth use 


1min F 


.09 Pg hfg 


gg(p f -p q ) 

L<W 2 


1 1/4 


(14) 


in equation (13) . 


A word about the scale parameter, L' , is in order at this 
point. In the form that we have given it, it has been named af- 
ter both Rayleigh who used it and Laplace who probably originated 
it. In this form it compares the size of the heater, L, with a 
fraction of either a capillary wavelength or a bubble departure 
diameter. In squared form it is exactly the Bond number. 


L' 2 


Bond number. 


g (Pf-Pq)L 3 
o L 


buoyant force 
surface tension force 


(15) 


By using L' instead of L' 3 in our studies we are laying heavier 
stress on its geometrical interpretation, than on its dynamical 
interpretation . 


If the influences of y and 0 C are not to be ignored then one 
more dimensionless group is required by the Pi theorem for each. 
e c is already dimensionless and N (recall equation (10)) is an 
appropriate group containing y . Thus we have 


Imax Imin 

Imaxp or Iminp 


f [L 1 , N, Pg/ (Pf — Pg) , 0 C ] 


(16 ) 


1 - 2 - 


We shall name, N, the induced convection buoyancy par ante ter 
since it compares gravity (or buoyancy) forces with the other 
forces in the following rather complex way: 


N = 


(inertia force) (surface tension force) 

2 1/2 
(viscous force) (buoyant force) ' 


3/2 


(17) 


The word "convection" in this context alludes to the flow induced 
by the viscous drag which is exerted by bubbles buoying upward 
from the heater. 

From a practical viewpoint we shall see shortly that it is 
advantageous to replace N in equation (15) with a scale parameter, 
i.e. , with a parameter that increases with L. It will also prove 
useful to replace it with a parameter that does not depend upon 
gravity. Accordingly we shall define an induced convection scale 
parameter , I , as 


I E-yW =yyfL CT / u 2 


(18) 


which characterizes the relevant forces in the following way: 
[ (inertia force ) (surf ace tens ion force ) ] 


I = 


viscous force 


(19) 


Thus we have as an alternative to equation (16) the follow- 
ing correlating equation: 


'Imax 


Plmin 


J-iLlCt-ft. -Kill _Li! p _» ■ . . -] 

qmax F <2min F “ f EL ' X ' Pg/(pf-Pg>' e c ] 


( 20 ) 


Some of the limiting behavior of this expression can be an- 
ticipated. We have already indicated that we do not expect to 
see Pg/(pf - Pg) exert a great influence, especially when it is 
small (or when the saturation temperature is much below the crit- 
ical temperature) . Likewise 0 C should normally have little 
importance with heaters that are not very small. 


Since L 1 and I are scale parameters, we anticipate that they 
will exert decreasing influence as heaters become large. The 
reason is that very large L 1 and I imply that buoyancy and inertia 
completely over-balance surface tension and viscosity. Once these 
forces have been over-balanced, further increases of L' and I do 
not correspond with any changes in the situation and we expect 
that 


Limit qmax ^min 

L)l->-<x. q ma x F ° r ^rninp 


constant 


( 21 ) 
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The low L' and I limits are harder to anticipate, and we shall 
defer a discussion of them to the last section. 

Equations (20) and (16) represent our basic correlation 
methods and we shall repeatedly apply them in the subsequent 
sections . 


14 


IV. FILM BOILING FROM HORIZONTAL CYLINDERS 


Theory 


The theoretical expression for the minimum film boiling 
heat flux on horizontal cylinders, equation (5), was developed 
in [10] on the basis of an assumed interface configuration shown 
in Fig. 6a. It was shown that the frequency of bubble departure 
was the following real number (iw) : 


(itjj) = 


kg 


P.F-P 


a - 


crk' 


crk 


p f +p g 


P f+ P g 2(p f +P g )R 2 


1/2 


( 22 ) 


where k is the wave number of the collapsing waves that feed the 
bubbles. It was also shown that the dominant wavelength, 
during film boiling could be obtained by maximizing this fre- 
quency. The result was: 


X d ~ 


2/3, 7T 
g( Pf~Pg ) 


H 


2R" 


(23) 


Finally, it was shown that by paraphrasing Zuber's [8] and 
Berenson's [17] arguments, equation (23) could be used as the 
basis for the prediction of q^-, equation (5) . Equation (23) 
was verified with a limited number of data for isopropanol and 
benzene in a fairly low range of R'. However, no attempt has yet 
been made to verify it under conditions of variable gravity. 
Equation (5) has been subjected to even less experimental verific- 
ation. 


The use of several transformations will put equations (22) , 
(23) , (5) into convenient dimensionless form. We first introduce 
a dimensionless wavelength: 

A E A d /A dF (24) 

where A d is Bellman's and Pennington's [21] wavelength for cap- 
illary F gravity waves in a plane, horizontal, liquid-vapor inter- 
face as might exist during film boiling above a flat plate. 


2/3 7T 

g(p^-p^) 


(25) 
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(a) ASYMMETRICAL INTERFACE WITHOUT WAVES 
ON THE BOTTOM 
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Fig. 6 Possible configurations of the liquid -vapor interface around 
horizontal cylindrical heaters during film boiling 


We next introduce a dimensionless frequency. 


ft = iwlcr/g 3 ( p f - p g ) ] - 1 / 4 


( 26 ) 


and we finally introduce the dimensionless radius, R' . Then equa- 
tion (22) becomes 


ft 


2 


1 

/3 A lw 


1 + 




p f p g 
Pf + Pg_ 


(22a) 


where in the present applications, (pf~pg) / (pf+p q ) - 1 in all 
cases. Equation (23) becomes y 


A LW R » 2^+ 1/2 (23a) 

where the subscript LW has been introduced to distinguish the 
Lienhard-Wong results from those of other authors we wish to con- 
sider. Equation (5) becomes, not equation (6), but 

^min/^minp = 0.515 f(R') (27) 

where the constant .515 arises because q m in-p is Berenson's result, 
equation (14), instead of Zuber's equation 12), and where 

f(R') = 


Data [10] in the range of dimensionless radius, R' , from .017 to 
0.25 showed that equation (23a) predicted wavelengths about 25 per- 
cent too low. In the discussion of [10] , the authors pointed out 
that this could probably be improved by including the minimum 
vapor blanket thickness, b, in the calculation. Two subsequent 
attempts were accordingly made to improve the accuracy of equa- 
tion (23a) . The first was by Siegel and Keshock [22] who assumed 
that the interface was shaped as shown in Fig. 6b 4 instead of Fig. 
6a and who retained b in their analysis. They obtained: 
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R ' 2 (2R 1 2 + 1) 


j./ ‘i 


(28) 


A SK 



R‘ 2 

+ (1 + b/R) -2 


(29) 


^This point might not be entirely clear in reference [22] which 
includes both the left-hand sketch in Fig. 6a along with the 
right-hand sketch in Fig. 6b. The analysis in [22] is based on 
Fig. 6b, however. 
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We can also modify equation (23a) so that it takes account 
of a finite vapor blanket. This is done replacing R with the 
radius (R+b) of the minimum interface radius in equation (23) 
and carrying out the nondimensionalization . The result is: 


A LW with 


(1 + -) 
R 



R 1 2 

+ 1/ [ 2 (1 + b/R) 



(30) 


While Siegel and Keshock did not evaluate b or make exper- 
imental comparisons, Baumeister and Hamill [23] succeeded in der- 
iving an expression for b, two years later. They used an assump- 
tion that the heat transfer must be maximum, subject to appropri- 
ate constraints, and their vapor "escapement" assumption differed 
from those-proposed in Fig. 6. Instead of taking the top of the 
interface to be wavelike they imagined spherical domes connected 
by annular passages. They obtained: 


A 


BH 


2 f\/T + 6R* 2 (1 + b/R) 2 - 1~ ) 
/3 it (_ Vr ,2 (1 + b/R) 2 


(31) 


Their analysis also permitted them to predict (1 + b/R) ? 


(1 + b/R) = exp{ 3 . 65 [kpAT/Rpgah^g] } 


(32) 


where hJl is a latent heat that includes a sensible heat correc- 
tion 

0.34c„ATl 2 
1 + — : h 


fg 


l fg 


fg 


(33) 


Comparison of Wavelength Results with Theory 

In an attempt to assess equations (23a) , (29 ) , (30) , and (31) 
for A, we made 233 still photographs during film boiling at over 
50 different conditions of gravity and radius [11] . Although we 
can scale lengths with an accuracy of about 3 percent, the wave- 
lengths exhibited wide variability. Most of them were obtained in 
acetone, but several. observations in methanol and isopropanol were 
included to insure that other fluids obeyed the same relation- 
ships . 


The heat flux was arbitrarily regulated to some value well 
below the peak and hopefully not far above the minimum. We found 
that, although Lienhard and Wong observed no effect of heat flux 
upon X in the range q/q m in - 4 , it was possible to distort the 
wave action and increase the merging of adjacent bubbles at higher 
heat fluxes. This is consistent with the studies Nishikawa et al « 
[24] who observed, not X, but the spacing of bubbles leaving 
wires in water over an eightfold variation of q up to the order 
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of magnitude of the peak heat flux, at R ' = 0.14. They found 
that bubble spacing increased about twofold in this range. 

Figure 7 shows six typical photographs of film boiling 
under a variety of conditions, arranged in order in increasing 
R' . Let us consider each of these: 

7a , Careful scrutiny of motion picture records in the R' 
range typified by this picture shows that the Taylor unstable 
wavelength briefly appears in the interface in the wake of depart 
ing bubbles . It is much shorter than the spacing between the 
bubbles on the wire and it very quickly merges into these bubbles 
Since the bubbles must grow quite large by successive mergers be- 
fore they have enough buoyancy to overcome the capillary forces 
that hold them to the wire, the bubble spacing is much larger 
than ^ 3 . We shall refer to this vapor removal process, which can 
be made out in Fig. 7a, as the "bubble-merger" mechanism. Grav- 
ity exerts little influence on the interface, and Fig. 6 b might 
possibly provide a more apt description of the interface than 
Fig. 6 a. 


7b . The Lienhard-Wong description fits almost perfectly 
here. Transverse surface tension, axial surface tension, and 
gravity are all in good balance. Vapor removal is accomplished 
by the orderly growth of Taylor unstable waves and their subse- 
quent collapse into departing bubbles. We shall subsequently 
refer to this as the "wave-collapse" mechanism of vapor removal. 

7c . The radius is smaller than in 7b but g is larger and 
the resultant R' is somewhat higher than in 7b. The heat flux 
is much higher than q m in anc 3 - r while a wave pattern is clear, the 
wavelength varies along the length of the wire. 

7d. The next value of R' is obtained by subjecting a small 
wire to a high g. The wave pattern is again clear and faithful 
to Fig. 6 a, but it exhibits almost 2 to 1 variability. 

7e . For a value of R' almost equal to unity, the relative 
diminution of transverse surface tension results in a deteriora- 
tion of the form of the wave pattern 

7f . The shaping effect of transverse surface tension is 
almost wholly lost in this case and estimates of wavelength were 
made mainly from the spacing of departing bubbles . 

The wavelength data are presented graphically in dimension- 
less forms in Figs. 8 and 9. Figure 8 also includes the low R' 
results from Lienhard's and Wong's stationary tests. The reader 
should remember as he scans Figs . 7 , 8 , and 9 , that while A in- 
creases monatonically with R' , X does not. Fig. 7d, for example. 
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DIMENSIONLESS WAVELENGTH, A 


1.5 


Experimental Data 

□ methanol 
a isopropanol 

v acetone 

0 benzene 

T variability of 

1 observation 


‘LW “ V R' 2 +l/2 


Open symbol* stationary test, I 
Solid symbol * centrifuge test, I <§i— 81 
•* reference [ 26 ] 


DIMENSIONLESS HEATER RADIUS, R' 

Fig. 8 Experimental wavelength on horizontal cylinders 
low R 1 range 




Dimensionless Wavelength , 


^ Eqn.(30)with (l+ p) for acetone 
\\\^\ Data from Fig. 8 

{ Centrifuge data for acetone 
with limits of variability 

A for Xi = 90% of 


Dimensionless Radius, R' 

Fig. 9 Experimental wavelength on horizontal cylinders — high R* range 
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corresponds with a large A, even though A is low at this high 
gravity . 

On the basis of the photographic results, we find that the 
Lienhard-Wong description is quite valid in the range 0.12 i R' 

< 0.9, but that it deteriorates in the range 0.07 < R' < 0.12. 

For R' 5 0.07, the bubble-merger mechanism is completely in con- 
trol of vapor dynamics . 

For R' > 0.9, the orderliness of the analytical model no 
longer exists in reality. Secondary wave motions and misalign- 
ment tend to disrupt the pattern. Fig. 6a appears to give only 
a reasonable first-order description of what is becoming an in- 
creasingly complicated phenomenon at larger R* . However, we were 
able to discern ragged waves up to R 1 > 2.5 and the Taylor theory 
still appears to account for the magnitude of these waves . 

We now wish to reconsider the three predictions of A : equa- 
tion (23a) or (30), equation (29), and equation (31). We have 
evaluated (1+b/R) on the computer for acetone, using equation (32). 
Fortunately this expression is only very weakly dependent upon AT 
in the ranges of interest. It was safe to assume AT - 333 e C 
since an error of 2 or 3 hundred degrees would not seriously 
change the predicted value of (1+b/R) . This computation is prob- 
ably not very precise, and it should only be viewed as indicative 
of the influence of b. 

Figure 10 shows A LW , A SK , and A BB together for comparison. 
The contribution of the vapor blanket thickness is slight in this 
case, as illustrated by the slight difference between the A B w 
curves including and neglecting b. The A SK and A BH curves both 
fall beneath the A LW curves and, as we can see in Fig. 9, the data 
fall on and above A LW . 

Thus, Ag K is probably based upon an erroneous description 
of the interface — at least in the range R' - 0.07. (A more de- 
tailed discussion of A in the range R' < 0.07 will be deferred 
until later in this Chapter.) The fact that the Baumeister- 
Hamill result is so low probably reflects the comparative crudity 
of their model for the configuration. They solved a much more 
complete overall heat transfer problem and had to make stronger 
simplifying assumptions. 

The broad variability of the wavelength data in Figs. 8 and 
9 is explainable if we keep sight of the basic assumption under- 
lying the prediction of both A B w and A SK . This is the supposi- 
tion that A is equal to A^, the wavelength for which the fre- 
quency, £2, of oscillation is a maximum. Figure 11 shows a plot 
of against A, based upon equation (22a) . This plot shows that 
the relationship is almost neutrally stable over a broad range 
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Dimensionless Radius , R 


Fig. 10 Comparison of predictions of A 







of A. The narrow range of frequency, 9 > 0.9fl max , subtends a wave- 
length range of 0.75A L ^ £ A <■ I.SOAlw c “" 

The implication is that there is a broad band of frequencies 
favored almost equally with ^ max . We should therefore expect to 
measure wavelengths in the range of, say, 0.75A LW to 1.6A LW . The 
choice of a 10 percent deviation of ^ below f2 max is of course ar- 
bitrary, but the use of this band in Fig. 9 appears to match the 
observed variability of data quite well. 


The Minimum Heat Flux 


The fact that we can predict X^ with reasonable accuracy 
under variable g conditions gives us added confidence in equation 
(6) for q m i n f which was developed in [10] on the basis of the 
prediction. We also know that we must distrust equation (6) be- 
low R'=0.07 because the mechanism upon which it is based is no 
longer valid in this range. 


In the present case we anticipate no influence of viscosity, 
there is no liquid-surface contact, and we expect that the major 
effects of elevated pressure have been accounted in the derivation 
of qininp • Thus equation (20) reduces to equation (8) for this 
case. Equation (6) is of this form, but Berenson's equation (14) 
should be used to express qminp instead of Zuber's equation (2). 
Thus equation (6) will change by a multiplying constant of 

[tt 2 1^73/60 ]/0. 09, and we obtain 


^min _ f 1.289 
q min F r ,2 (2R' 2 +1) 


(34) 


Equation (34) can also be redeveloped to incorporate the va- 
por blanket thickness introduced through X^. To do this we must 
redo the analysis in reference [10] . This is straight forward and 
the result is not surprising. It is ( cf . equation (29) ) 


qmin 

^minp 


1.289 

R ,2 [2R ,2 +(1+ ^.) -2 ] 
R 


1/4 


(34a) 


On the basis of computations of (1+b/R) for typical situations, 
using equation (32) , we can conclude that the blanket thickness 
will not alter qmin/^minF ky more than a few percent in the worst 
cases. Therefore, we shall only consider predictions in the form 
of equation (34) and hereafter ignore equation (34a) . 


Before attempting to compare equation (34) with data, it is' 
well to consider a serious experimental difficulty in obtaining 
q m in for horizontal cylinders. References [10], [2], and [11] 

refer to the difficulties of eliminating end effects which cause * 
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the f ilm-to-nucleate boiling transition to occur at heat fluxes 
above the true q m in* Very recently Kovalev [25] devised an ex- 
perimental trick for determining the true qmin* He bent the ends 
of his horizontal cylindrical heater 90° into the vertical posi- 
tion and led them out of the bath. In doing so he removed the 
end effect problem and obtained qmin data that are probably re- 
liable. His results are lower than other q m in data that have 
been reported at comparable conditions . 


Another recent study by Grigull and Abadzic [26] gives q m j_ n 
data for CO 2 and Freon-13 at extremely high reduced pressures, 
which we shall wish to consider. However, there is generally a 
paucity of reliable q m in data for cylinders , owing to the prob- 
lems of end effects and surface impurities that affect so many 
sets of data. In particular, there is a lack of usable data in 
the range R ' > 1.0. 

Data from [10] , [25] , and [26] are compared with equation 
(34) in Fig. 12. The recent data of Kovalev and Grigull-Abadzic 
are lower than those from [10] — probably owing to even greater 
care taken to avoid end effects. The highest R' points in both 
the CO 2 and Freon-13 series were measured so near the critical 
temperature as to reflect some of the anomalous behavior assoc- 
iated with the critical point. Therefore we have entered these 
points with dashed lines . The Grigull-Abadzic paper also in- 
cluded one photo made at low R' and at 96 percent of the critical 
pressure, from which we extracted an approximate data point for 
inclusion in Fig. 8. 


Since the multiplying constant in equation (34) is empirical 
it should probably be chosen lower than it was in [10] to give a 
better representation of the ultimate hydrodynamic minimum film 
boiling heat flux. Accordingly we recommend 


Qmin 

qminp 


0.0217 

R , 2 (2R» 2 +i) 


1/4 


(35) 


for the real minimum as suggested by the data in Fig. 12. Equa- 
tion (35) is included in Fig. 12. 


From a practical viewpoint, however, equation (35) will al- 
most always predict lower qmin' 3 than will be observed, since 
every normal experimental difficulty will tend to give values 
greater than the ultimate hydrodynamic minimum. For this reason 
Berenson [17] spoke of a minimum, minimum heat flux. Indeed, 
Kovalev's experiments cast the same kind of ugly light on q m i n 
measurements that the early surface roughness measurements of 
Gorty and Faust [27] cast upon nucleate boiling. In either case 
the variable of interest appears to become swamped by capricious 
system variables . 


/ 
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dimensionless minimum heat flux, q m j n / 


04 

O 



Fig. 12 Predicted and observed minimum heat fluxes on horizontal cylinders 




The Limits of Large and Small R 1 


Large R* . The clean wave-collapse mechanism described in 
[10] and displayed in Figs. 7b through 7d or 7e loses its coher- 
ence as R' exceeds 0.9. Nevertheless, the observed "wavelengths" 
still obeyed equation (23) up to the limit of our experiments at 
R'=2.52 — at least in-so-far as wavelengths could be identified. 

In this range, the measured values and equation (23) both approach 
Xdjp as a limit. 

There is some evidence that A< 3 F is approximately the correct 
limit for even larger R' . Pomerantz [28] measured heat transfer 
coefficients during the film boiling of Freon-113 on 3/16 in. dia- 
meter tubes, under variable gravity. He also gave three ranges of 
bubble spacings (which probably approximate Taylor wavelengths at 
his large R ,l s and relatively low heat fluxes). They reduce to 

A(R' = 2.27) = 0.5 to 2.32 
A (R ' = 4.55) = 0.59 to 1.0 
A(R' = 5.9) = 0.65 to 1.18 

All three of these values scatter widely about A = 1, or A<j = A^p. 
He also claimed an average value of A = 1.1 for his bubble spac- 
ings, which is consistent with the trend of our wavelength data. 
The early photographic study of Westwater and Santangelo [29] pro- 
vides one very clear photograph of methanol in film boiling on a 
horizontal tube from which we obtain 

A(R' = 3.03) = 0.80 to 1.1 

This also corroborates the trend at large R' . 

These observed wavelengths at large values of R* are all 
probably low by a small but unknown amount owing to foreshorten- 
ing, since the waves no longer align in the axial plane. Unfor- 
tunately, no one has yet documented the way in which three- 
dimensionality develops at large R' . 

The minimum heat flux prediction — equation (35) — will be 
valid until these three-dimensional effects start to be signifi- 
cant. This will probably occur above R' =2. Once the wave 
pattern ceases to align on the axis, more bubble releases can be 
accomodated on top of the wire, and the minimum heat flux should 
exceed the prediction. We have already noted that qmin/^minp 
should level out to a constant value (recall equation (21)). In 
the subsequent sections of this report we shall show that this is 
clearly the case for qmaxAlmaxp ^- n several configurations and for 
‘Imin/qminp in at least one other case. 
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But. there are no reliable q m i n < ^ a ‘ t:a available fox large 
radius horizontal cylinders, and the limiting behavior is simply 
not known at this time. A recent and very extensive survey of 
film and transition boiling by Jordan [30] , for example, provides 
only one additional q m i n data point that we might have used. 

This was quoted from [29] and was far too high to have represented 
a true minimum. 

Small R 1 . The transition in the range 0.07 < R" £ 0.12, 
between the "bubble-merger" and the "wave-collapse" mechanisms is 
much more severe and far more distinct than the gradual deterior- 
ation of the wave-collapse mechanism at large R' . The following 
mechanism can be established for this transition: 

A study of many photographs shows that bubble departure 
takes place roughly as shown in Fig. 13. The radius, X^/16 , 
shown in the sketch appears to be characteristic of this config- 
uration. 

This sharp corner can exist as long as it is not opposed by 
an equally strong counter-curvature of the interface around the 
cylinder. But for small wires the counter-curvature will prevent 
the corner from indenting inward and nipping off the bubble. At 
this point the departing bubble will have to spread out before it 
can depart, and merge with its neighbors, if R is sufficiently 
small. 

It is now possible to compare A^/16 with the radius of the 
wire at the observed transition from the wave-collapse mechanism 
to the bubble-merger mechanism. We have the condition: 

R' 2 0.07 for wave-collapse mechanism (36) 

which becomes „ 

R > 0^07 _ 0 . 00644 

A d F 2/3ir 

But from equation (23a) 

x d 
X d F 

R =0.07 

Therefore condition (36) takes the form 

R „ 0.00644 _ 1 

X d - 0.0986 “ 15.3 

which corresponds very well with photographic observations . 


0.07 

-2 ^ 


= 0.0986 


+ 1/2 
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The problem of developing a description of vapor removal 
from small cylindrical, heaters in film boiling at low gravity 
remains to be done. So too does the prediction of q m j„. When 
they are accomplished such analyses will focus upon a force bal- 
ance between surface tension and buoyancy. They will not be 
based upon considerations of Zuber's capillary waves. 


Conclusions 


1. The best prediction of A that we presently can make is: 


A = nr 2 — y (30) 

1 / R 1 ^ + 1/ [2(1 + b/R) 2 ] 

where (1+b/R) can best be estimated at present using equation 
(32) . The variability of actual data will be slightly less 
than +60 percent and -25 percent because the wave motion ap- 
proaches neutral stability over this range. 

2. It is not relevant whether A is changed by varying, R, g, or 
the liquid being boiled. The single parameter, R' , correctly 
accounts for all of these effects upon A. 

3. Our best prediction of the absolute minimum q m ^ n for horizon- 
tal cylinders is: 

^min/^max 

4. For— R-'_-< 0.07 no mechanism has yet been formulated to describe 
the vapor removal process or q m j_ n . 

5. For R' >2 the hydrodynamic mechanism for the vapor removal 
process is unclear and q m ^ n is not known. Reliable q m ^ n data 
and extensive photographs are needed in this range. 


0.0217 


R' 2 (2R' 2 + 


1 ) 




(35) 
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ON HORIZONTAL CYLINDERS 


V. 


^max 


Prefatory Remarks 

Without doubt the most thoroughly studied pool boiling con- 
figuration todate has been the horizontal cylinder, and rightly 
so, since so many cylindrical heaters are used in practice. The 
quest for the influence of gravity on q max ^ as more often than 
not focused on this geometry (see, e.g., the review by Siegel 
[31] ) , and a great deal of data exist for cylinders under a var- 
iety of conditions . 

To some 4 or 5 hundred existing q max data we have added 
about 440 of our own. These data provide an extremely complete 
verification of our method of correlation and of the prediction 
that we shall offer. We shall first develop the prediction and 
then consider the available data. 

The prediction will be based on the three assumptio.ns sug- 
gested after equation (9). Of these, the second two assumptions — 
that 0 C and Pg/(Pf~Pg) are unimportant — have been adequately dis- 
cussed. The assumption that y can be neglected will be largely 
justified ex post facto . It is not unreasonable because the 
flows induced by bubbles rising above a cylinder will not be 
drawn into the bubble departure pattern. The next Chapter will 
treat a different situation where this is not the case and y is 
most important. 

Under these assumptions we can anticipate that the general 
functional expression for qmax r equation (20), will reduce to 
the simple form 


■ - = f (R' ) (8) 

<3max F 

This is the form that our predictions will have to take. 

The previous Chapter showed that we should be alert for 
three regimes of hydrodynamic behavior. In the first R' << 1 
and capillary forces completely overwhelm buoyant forces. For 
the moment this regime will be considered to lie outside of our 
range of interest since it cannot be treated by any adaptation 
of the hydrodynamic theory. 

The middle regime is one in which capillary and buoyancy 
forces exert comparable influence and their ratio (R')2 should 
be on the order of unity. This is the range in which we are 
presently interested, and any theory we reach should reflect a 
distinct interaction of these forces . 
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In the third regime — R' >> 1 — the influence of capillary 
should cease to have any role in shaping and scaling the process . ■ 
It is in this regime that theory should yield a result in the 
form of equation (21) . 

The primary assumption that must be made in formulating the 
problem is the idealization of the vapor removal configuration. 

Our vapor removal hypothesis is sketched in Fig. 14. We shall 
consider a horizontal cylindrical heater submerged in saturated 
liquid. Vapor is generated around the cylinder and rises in the 
direction of acceleration of the system. We imagine that there 
is a plane above the cylinder, where the bubbles from the upper 
half and lower half of the cylinder merge and form vapor columns . 
On this plane, the liquid is supported by the vapor in a Taylor- 
instable interface. The jets will be spaced on the Taylor- 
instable wavelength, A d F- This wavelength is chosen in prefer- 
ence to the wavelength for a cylindrical interface since the con- 
figurations shown in Fig. 6 were not replicated in our observa- 
tions of boiling at the peak heat flux transition. 

Figure 15 includes five typical photographs of boiling just 
below the peak heat flux transition, over the range 0.0695 1 R'i 
2.68. These are typical of some 640 such pictures that we made 
to serve as a basis for the several suppositions made thus far 
about the character of the process, and for those that will be 
made subsequently. 


q m ax f° r on fh* 2 Order of Unity 


We shall base the analysis upon consideration of a "unit 
cell" of one wavelength, A d p. If Ajj is the area of the heater 
surface within a cell and Ag is the cross-sectional area of a 
vapor jet, then 

Ag it (R+<$ ) 2 _ (R+6) 2 

A H ~ 2-n-RA d 2RA d ' 

The heat flux, q, from the heater during saturated boiling is bal- 
anced by the latent heat carried away in the vapor bubbles . Since 
all the bubbles merge and form vapor jets with area, Ag, we can 
write an energy balance in the form 

q = p g hfg Ug -2 = P g hfg U g ^RA^ 

where Ug is the velocity of the vapor in the jet. The heat flux 
reaches qmax when the velocity of the vapor jets reaches the 
critical value, Ugc, at which they become Helmholtz- instab le . 
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H-2(R+8) -H 



(b) Cross-sectional view of a vapor jet. 


Fig. 14 Model for the peak pool boiling heat flux on horizontal 
cylindrical heater, (R 1 on the order of unity) . 
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Fig. 15 Photographs of nucleate boiling at q just below q max 








Therefore, 


Pg kfg Ug C 


(R+6) 2 
2RX d 


(37a) 


Next, we must determine Ug C . We shall assume, as Zuber did, 
that the Helmholtz-instable collapse of the jet occurs as a result 
of a major disturbance arising from Rayleigh's capillary instabil- 
ity. According to Rayleigh [32, p.473], the minimum wavelength 
that will be capillary- instable is one circumference, or 2tr(R+6), 
in length. This will be the major disturbance in the jet and it 
will be Helmholtz-instable when Ug reaches a value Ug C such that 
[32, p . 462] 


2 it (R+S) = 2tt 


Pf + Pg g 

P f Pg _f U f +U gc^ 


Thus, if we neglect the velocity, Uf, of liquid in-flow, and em- 
ploy the equivalent assumption, Pg«Pf, 


f Pg(R+S) 

Substituting equation (39) into equation (37) gives the 
equation for the peak heat flux on a horizontal cylinder 


Plmax ~ ~ 


(R+8 ) 2 

2Rl d -p 


Pg h fg 


P g (R+6) 


While the radius of the cylinder and the physical properties are 
known, the vapor blanket thickness, 6, remains to foe determined. 
If we define a dimensionless vapor blanket thickness 

A 5 <$[g (pf-p g )/cr] 1 / 2 (41] 

and combine equations (25) and (41) with equation (40) , we get 


%ax 


"JThfg [ag(p f -p g )] 1/4 ]|^-= 


(R’+A) 


The first, group on the right hand side of equation (42) contains 
a constant which has been balanced against the second term so 
that it matches Zuber' s flat plate prediction, q^axp • Thus 

Ptmax _ 6 (R-M- A ) 3 ( 43 ' 

qmax F “ 1T 2 /3 - RI 

where A remains to be determined as a function of R' , with the 
help of experiments. 
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Prediction of q max for R' >> 1 


We have assumed that vapor jets with a radius of (R+<5) are 
spaced on the wavelength, . But A^p decreases with increasing 
gravity, g, so the spacing between jets will decrease as R' in- 
creases. According to the previous analytical model, the jets 
will touch one another when (R'+A) reaches 2 tt/ 3/2 . However, we 
did not see this happen in the present investigation of pool boil- 


ing, 


Therefore, the preceding model should fail at some (R'+A) 


smaller than tt/3. An illustrative sequence of sketches showing 
the supposed variation of the jets with increasing R 1 is given 
in Fig. 16 . 


When R' is small, the unstable wavelength, A 3 , is greater 
than twice the jet diameter, 4(R+5), as shown in Fig. 16a. As R' 
increases, the difference between A 3 and 4(R+S) becomes smaller. 
We believe that, when R' is very large, the open space between 
two jets will be equal to 2(R+<5) and the jet configuration will 
remain unchanged. The jets will stay as shown in Fig. 16c in- 
stead of moving any closer as indicated in Fig. 16b. This is 
analogous to Zuber's assumption that the jet diameter on a flat 
plate is equal to the open space between two jets, and our photo- 
graphic studies verified that it is a realistic assumption. 


When R' becomes very large, disturbances with the wavelength 
A 3 , will exist in the horizontal vapor-liquid interface, but they 
will be much smaller than the jet diameter. The vapor jets will 
pick up these disturbances from the intervening liquid-vapor inter 
face. It is probably this disturbance that triggers Helmholtz- 
instability instead of the longer wavelength of the Rayleigh- 
instability. Thus we shall assume that the jet diameter is 
2(R+6), but that the distance between centerlines of jets is 
4(R+<5) instead of A 3 F . Finally, we use A 3 „ instead of 2 tt(R+ 6) 
for the left-hand side of equation (38) ana obtain 


Ugc 



(44) 


in the place of equation (39) . 


Combining equations (25) and (44) with equation (37) gives 


^max 


tt , 4 /crg(p f - Pcf ) 

24 Pg f *V P 2 

' g 


3 3 / 4 (R'+A) 
ir R* 


(45) 


5 

It was observed by Vliet and Leppert [33] that, in a strong 
forced convection boiling system, no jetting occurred. The vapor 
appeared instead to cover the whole wake region of the cylinder. 
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or 


%ax _ 3 3 / 4 (H 1 + A ) 
^maxp it R ' 


( 46 ) 


Equations (43) and (46) should each give the same peak heat 
flux at the point where one configuration shifts to another. Thus, 
we can write, for this particular value of R' , 

3 3 / 4 (R'+A) = 6 (R'+A) 2 / 3 (4?) 

ir R' tt2/3’ R' 

This gives, for the transition between the two analytical models, 

R' + A = 4.28 (48) 


Finally, the expression for q ma x at ver Y large R' should 

become 

qmax = ^±4 = constant, R' + A * 4.28 (46a) 

'Jmaxp it R ' 

where the dimensionless vapor blanket thickness still must be de- 
termined. Once A is known at (R'+A) =4. 28, the constant in equa- 
tion (46a) can be fixed. 

The jet diameter reaches 1^/2 at (R'+A)=ir/3/2 according to 
the model for R' on the order of unity. Beyond this size the 
open space between the two jets must become less than the jet 
diameter (or, conversely, the jet spacing might start to increase 
beyond X(j) until (R'+A) reaches 4.28. Thus (R'+A)=4.28 is act- 
ually the upper bound for a transition which might occur gradually 
over the range tt/3/2< (R' +A) <4 . 28 . 

The two predictions — equations (43) and (46) — require meas- 
urements of 6 before they can be completed, and they require con- 
siderable qmax data for verification. But we shall first consider 
the correlation of qmax data in accordance with equation (8) . 


Correlation of q m ax Data 


A total of about 900 qmax data were collected, both from 
the present experiments and from other sources . These observed 
q m ax were divided by Zuber ' s qmaxF ' which was computed with some 
care in [13] and is plotted for typical fluids in Fig. 17. 

Figure 18 displays the present data, as well as data of 
Cumo, Farello and Pinchera [34] and Lienhard and Watanabe [2] , 
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Fig. 18 Comparison of prediction with experimental data 
R 1 < 0.50 


for 



in the range R'<;0.5. It includes data for acetone, benzene, car- 
bon tetrachloride, isopropanol, methanol, and water over wide 
ranges of pressure, gravity, and heater size. The variability of- 
the data is on the order of ±20 percent except when R'cO.15. 


The correlation curve from reference [2] has also been in- 
cluded as a dashed line. This curve was originally presented on 
different coordinates, and in placing it in Fig. 18 we have fav- 
ored the original data - for acetone, slightly. The small irreg- 
ularities in this correlation are probably the result of minor 
changes in the transition mechanism, too small to be accounted 
in our relatively coarse analysis. 


Figure 19 shows our data for a larger range of R' . Data 
given by Adams [35] , Carne [36] , Costello and Heath [37] , Cumo, 
Fatello and Pinchera [34] , Frea and Costello [38] , and Pramuk 
and Westwater [39] are also included. All these data correlate 
consistently in gmax/^maxp versus R' coordinates. 


Figure 20 shows our data for high R' and additional data 
from references [35] and [37] . The data appear to approach a 
constant value of grnax/^max'F ^ or ^ ar 5 e as we anticipated. 

This tendency was also observed by Lienhard and Keeling [15] in 
their study of horizontal ribbon heaters. Lyon's [18] data for 
boiling liquid N2 and 02 are also included in Fig. 20. They all 
fall within the scatter of ±20 percent. 


Figure 21 is a general view of the results of the correla- 
tion. A variety of very low R' data obtained by Siegel and How- 
ell [5] under reduced-gravity conditions are also included. The 
correlation is shown to be very successful with the existing data 
for all R’>0.15, but when R'<0.15, the data show a wide scatter. 


The success of equation (8) in bringing together so many 
diverse data indicates that the three assumptions upon which it 
is based must be correct. In particular, it is most heartening 
to discover that neither surface condition nor viscosity has 
proved to be important — at least for R'>0.15. 

In the present study, we have collected data for many 
smooth and very clean surfaces, in which only 0 C might have been 
subject to variation. Only in the very low R' range do we dis- 
cover a failure of the simple correlation scheme. This might 
imply either technical difficulties in obtaining the experimental 
data or that equation (8) is incomplete in the low R' range. In 
particular 9 C or I might be important independent variables. 

The photographic studies (see, e.g.. Fig. 15a) showed that bubble 
growth and merger mechanisms govern the vapor removal process at 
very small R' . In this range, both viscosity and contact angle 
may become very significant. 
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Experimental Determination of A(R f ) and 
Completion of Predictions 

We have now verified the functional equation for q max by 
showing that data correlate in the dimensionless coordinates of 
this expression. Our prediction equations are also of this form 
so we can now be fairly optimistic about their success . 

To complete the predictions, however, it is necessary to 
evaluate A=A(R'). This was done by measuring the vapor blanket 
thickness from the photographs. Since the blanket has to be 
viewed from the side, § cannot be measured on the horizontal di- 
ametral plane. Hence, we had to approximate 6 with a measurement 
of the blanket thickness on the bottom part of the cylinder as 
shown in Fig. 14b. 

Figure 22 displays the measurements of 6. In an attempt to 
compensate the fact that the blanket might not be fully developed 
at the bottom, we have reported the range from the approximate 
mean, to the maximum size of the bubbles comprising the blanket. 

A data point measured from an early photograph by Westwater and 
Santangelo [29] is also included. The uncertainty of these meas 
urements is on the order of ±20 percent. No attempt was made to 
measure 6 in the range of R'<0.15. 

If we had to combine the envelope of the data from Fig. 21 
with equation (43), we can specify a range of A, consistent with 
the qmax theory, for each value of R* . The result is plotted in 
Fig. 22. Virtually all the measurements of 6 fall within this 
envelope . 

The "best fit" empirical equation, based on both this en- 
velope ana the data for 6 , is 

A = [2.54 R' + 6.48 R' exp(-3.44 /R r ) ] 2 / 3 - R' (49) 

in the range: 0.15<R'<3.5. The form of - equation (49) is chosen 
so as to give a neat form for the qmax/Qniax-c' equation that we 
shall develop next. F 

The substitution of equation (49) in equation (43) yields 


Qmax 

Qmaxp 


0.89 + 2.27 exp (-3.44 /R') 


(50) 


Equation (50) is plotted in Figs. 18, 19, 20, and 21 and it shows 
very good agreement with the existing data. 

The upper limit of equation (50) — the transition point be- 
tween the two models — is determined by equation (48) . With the 
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Fig, 22 Experimental determination of the dimensional vapor blanket thickness 
as a function of R'. 








help of equation (49) , we obtain 
R' = 3.47 at 


R' + A = 4.28 


be fixed 


The constant value for q m ^x/qmax]? e Q ua tion (46a) can now 
ed with the help of equations (51) and (49) : 


Qmax 

Qmaxp 


3 3 /4 4.28 
ir 3.47 


= .894; 


R' > 3.47 


Thus the proper upper bound on equation (49) is R'=3.47, and the 
relation between A and R' , for R'>3.47, can be found by combining 
equation (46a) with equation (52) to get 


A = 0.233 R’ 


R' > 3.47 


Equation (53), plotted on Fig. 22, is consistent with the few 
measured values of the vapor blanket thickness in this range. 

Figure 23 shows equation (43) for various values of A and 
the semi-theoretical equation for qmax/Qmaxp * Actually equation 
(50) dif fers negligibly from equation (52) * for R'>3 . 47 . Therefore 
equation (50) can properly be used for the entire range of R'>0.15. 

We have presented two models for the q ma x transition. Act- 
ually, there are small aspects of the transition mechanisms which 
are ignored in these relatively simple descriptions . The range 
of R' between 1.72, at which (R 1 +A) =tt/3/2 , and 3.47 is, for ex- 
ample, a transition range between two models. As R' increases in 
this range, the jet-spacing becomes decreasingly dependent on 


Conclusions 

1. A hydrodynamic model is developed for the prediction of q m ax 
on horizontal cylinders. The resulting expressions are 


Qmax 

Qmaxp 


6 , ( R ’+A . ) . 3/2 ; r. + A £ 4.28 


VT R' 


qmax 

Qmax-r 


o 3/ 4 TJ | i A 

- = constant; R' + A > 4.2 8 

7T R 1 


2. The dimensionless vapor blanket thickness. A, is determined 
by an empirical equation 

A = [2.54 R' + 6 . 48 R' exp(-3.44 /r^ - ) ] 2 / 3 -r' ; R'<3.47 
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Fig. 23 Theoretical dimensionless peak heat fluxes for various values of A, compared 
with final prediction. 


u> 







from which we find that R'=3.47 at R'+A=4.28. Therefore, 

A = 0.233 R» ; R* > 3.47 

3. About 900 q max data for cylindrical heaters of various sizes, 
under different pressure and gravity conditions, and in a 
variety of liquids were correlated successfully on gmax/^maxp 
versus R' coordinates. 

4. The resulting equation for q max on horizontal cylinders is 

‘Smax/^maxp ~ 0.890 + 2.27 exp (-3. 44 /R 1 ) ; 0.15 < R' 

5. When R'<0.15, the q max data start deviating from the predic- 
tion and scatter widely on the qmax/^maxp versus R' coordi- 
nates. When R'<0.07, the hydrodynamic instability theory of 
transition is no longer applicable. If a transition of any 
kind still exists in this range, it is one in which capillary 
forces are dominant . 
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VI. q ma x 0N HORIZONTAL RIBBONS 


Th e R o l e' o f N o r I 


Viscosity appears to exert no influence whatever on the peak 
heat flux on horizontal cylinders. But there is evidence that it 
might be very significant in certain more susceptible configura- 
tions. Borishanski [16], as we observed in Chapter III, found 
some combined influence of viscosity, buoyancy, and surface ten- 
sion upon q max . 

A much stronger indication of such influence was given by 
Costello, Bock, and Nichols [40] in 1965. They observed the burn- 
out of ribbons mounted in water as shown in Fig. 24. The rising 
bubbles induced a transverse flow which altered q m ax* The smaller 


Fig. 24 



Configuration studied by Costello, Bock and Nichols 


ribbons burned out at higher heat fluxes than the larger ones , and 
the lowest heat fluxes were observed when vertical sidewalls were 
placed along the edges of the ribbon to block the liquid inflow. 

Costello et al . also measured rates of liquid inflow, and 
they suggested that a superposition of forced-convection and pool- 
boiling effects might account for heat transfer in this case. 
However, they did not provide a correlative procedure and they 
gave data for only three situations. 

We have repeated their experiment in a slightly modified 
.configuration over a wide range of sizes, gravities, and pressures, 
in five fluids. We used the heater on the right hand side of 
Fig. 3 in the test capsule on the left. This yielded the cross- 
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sectional configuration shown in Fig- 25. 

In this configuration we would expect the following modif- 
ications in the basic correlation equation, equation (16) or (20) . 
The dimensionless width, W’ , should replace L' , where W’ is simp- 
ly an L' based on L=W. Since we are using common nichrome rib- 
bons in a range of fairly substantial W' , the effect of contact 
angle should be small and we shall ignore it. Our pressure range 
is low so the influence of Pg/(pf~Pg) should also be negligible. 

Finally, a secondary parameter is introduced into the prob- 
lem by the limited_ width , W c , of the capsule. We shall assume 
that, since recirculation results from the confinement of the 
capsule width, W c should appear in a second induced convection 
scale parameter, I c . Since the capsule is a good deal larger 
than the ribbon, recirculation should not exert a large influence 
on Ptmax/ but it cannot be ignored. Accordingly, we can expect 
to correlate data with 


^max 

^maxp 


f(W', I, I c ) 


(54) 


If we wish to use equation (20), which incorporates Bori- 
shanski 1 s induced convection buoyancy parameter, N, then it is 
not possible to introduce Wc in a second N, since N is not a 
scale parameter. We must instead use an aspect ratio, W/Wc, so 


g 51 * " = f (W', N, W/W c ) (55) 

%ax F 

The primary objective of the work described in this chapter 
is to learn whether or not any induced convection influence act- 
ually becomes significant. We do not anticipate a strong influ- 
ence of I c but its effect will also be investigated. If we are 
wrong in looking for an influence of I or N then equations (54) 
and (55) will reduce once again to equation (8). 


Development of Correlation Surfaces from Data 


Carefully cleaned nichrome ribbons were cut from 0.228, 
0.051, or 0.025 mm stock and stretched out on the mounting block. 
The capsule was charged with reagent grade acetone, methanol, 
benzene, isopropanol, or double distilled water, and mounted on 
the centrifuge. During a run, the liquid level was held in the 
range 2.5 to 3.8 cm above the ribbon surface and noted to make a 
hydrostatic head correction, and the temperature held to within 
a degree of saturation. The angular speed of the centrifuge 
was read, and (with the preheater turned off) the cap- 
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sule pressure was recorded. Finally, the power supplied to the 
ribbon was increased until the peak heat-flux transition was ob- 
served, and there it was read. 

About 874 observations were made — roughly 419 in acetone, 

262 in methanol, 104 in isopropanol, 81 in benzene, and 8 in water. 
These raw data are tabulated fully in reference [14] . The prob- 
able errors of the variables, computed in [14], are about ±6 per- 
cent for q ma x' ±4 percent^for W' , and ±2 percent for I. The var- 
iability of observed q maK values was on the order of ±15 percent, 
as is typical of q max data. Data were measured over a reduced 
pressure range from 0.0016 to 0.0246, and a range of gravity from 
1 to 87 times earth normal gravity. 

All of these data and Costello's three points were plotted 
twice: first on q max /qmaxp versus W' coordinates for comparative- 

ly narrow ranges of I; then on qinax/^maxp versus I coordinates for 
ranges of W' . Fig. 26 shows a typical example of one of these 
crossplots . ° 


Figure 26 reveals some things that were generally true of 
all the data. The great majority of the data for any substance 
clustered within ±15 percent of a mean surface through them, and 
a geometrically similar family of surfaces could be drawn through 
all of the substances. Methanol was most representative of the 
substances used. Acetone and water presented the highest qmax/ 
qmaxp values — about 25 percent above methanol — and isopropanol 
the lowest — about 30 percent below the methanol. The benzene data 
lie between the methanol and acetone data. 

This is consistent with the fact that we have failed thus 
far to account for the remaining weak parameter I c . We shall re- 
turn to this point briefly, but let us continue for a moment as 
though the influence of I c could be neglected. 

The correlating surfaces that result from the crossplots are 
each presented in sets of contours. The correlation function 


/r 

The method of correlation used in reference [14] was based on a 
Law of Corresponding States correlation instead of equation (54) . 
This resulted in a doubling of those errors introduced by the Law 
of Corresponding States and an awkward scaling of q max . The pres- 
ent correlation overcomes these difficulties but required replot- 
ting the data instead of using the curves in [14] . The effective- 
ness of equation (54) was, of course, made clear in Chapter V. 

The shape of our resulting surface differs somewhat from that 
plotted in reference [14] . 
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Dimensionless 





f (W 1 , I) is presented in Figs. 27 and 28 which give q max /q ma xp 
versus W' and qmax/^maxp versus I contours, respectively, as* ob- 
tained from the 21 crossplots. This surface generally follows 
the mean of the data for methanol and is consistent with trends 
that can be identified individually in each of the remaining sub- 
stances . 

The correlation in terms of N was obtained by transforming 
th e cu rves given in Figs. 27 and 28 with the help of equation (18), 
I=/NW' . It is given in a single plot of qmaxAlmaxp versus W' , 
with N as parameter, in Fig. 29. 

Figures 27, 28, and 29 are- limited in their applicability 
to just methanol. More specifically they are restricted to the 
present configuration in a narrow band of I c . Table 1 gives I c 
ranges for the fluids used here. 


Table 1 I c for the Present Tests 


Liquid 

I c for atmos- 
pheric pressure 

I c for minimum 
pressure 

Represent 
ative I,, 

isopropanol 

2150 

948 

2000 

methanol 

3385 

2375 

3200 

benzene 

4100 

2190 

3800 

acetone 

4650 

3300 

4500 

water 

7520 



7520 


Figure 30 shows constant W' curves for acetone and isopropanol 
to show how I c influences the surface. This figure is in the 
same form as Fig. 28, and two lines in it have been borrowed 
from Fig. 28 to illustrate what the relative influence of I c is. 


Discussion 


The correlation surfaces reveal the strong influence of 
the parameters I and N, that we anticipated, and they also show 
that this influence vanishes when the scale parameters, I and W' , 
become large. This was to be anticipated in the context of equa- 
tion (21) , since buoyancy completely overbalances both capillary 
and viscous forces as the scale is increased. Thus we find for 
the methanol data that: 
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Fig. 27 q m _„/q. 


max' ^max t 


vs. W* contours for 9 values of I 
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limit ^ max - 1/2, large I and W f (56) 

Imaxp 

This limit appears to be valid for all I>1500, and W' on the 
order of 50, depending upon the value of I. Since N is not a 
scale parameter we cannot propose a proper criterion in terms of 
it. However, the various lines of constant N seek qniax/ < 3maxT? = ^/ 2 
as an asymptote in Fig. 29. 

That this limiting q max is less than the flat plate value 
might reflect either or both of two factors: The primary factor 

is the finite size of the centrifuge capsule which doubtless re- 
sults in a limiting high Reynolds number circulation. This circ- 
ulation in turn affects the peak heat flux. A secondary factor 
might lie in the fact that Zuber's equation has not really been 
subjected to broad testing in a proper flat plate configuration, 
and some error in the constant, ir/24 , could conceivably contri- 
bute to a deviation of the present data from the supposed flat 
plate equation. 

When W' falls below unity, surface tension assumes very 
strong control over the peak heat flux transition and q max rises 
sharply. We saw in Chapter V that this was also true as 2R' — 
the comparable scale parameter for cylinders — fell below 0.40. 

Figure 30 shows that while the induced convection effect 
is pronounced, it is also susceptable to the geometry of the 
container as reflected in I c . The surfaces for acetone fall 25 
or 30 percent above those for methanol while those for isopro- 
panol fall about 30 percent below. 

Throughout the preceding discussion we have made no mention 
of the thinness of the ribbons used in the tests , however the 
method of correlation involved no consideration whatsoever of 
heat conduction or capacitance effects. These effects can influ- 
ence q max as the ribbon becomes thin and they should be discussed 
further . 

The low thermal capacity of thin ribbons led to surface 
temperature oscillations, owing to the alternating current power 
supply. This required that we change to direct current supply 
when the ribbon thickness became small. The criterion for this 
change was based on the heat conduction solution for this prob- 
lem given by Switzer and Lienhard [41] . Overlapping ac and dc 
experiments verified that we did not suffer difference in results 
at the point of transition. 

. It has long been observed that very thin ribbons can actu- 

ally melt out at heat fluxes less than the hydrodynamic limit. 
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q , if they lack the protection of enough "thermal mass." 

This burnout mechanism is discussed at length in [42] and a cri- 
terion is given for determining the limiting thickness. We took 
care to keep our observations within this limit. 


Conclusions 


1. Induced convection can exert a strong influence on the peak 
pool-boiling heat flux, if the configuration is one which is 
susceptible to it. 

2 . The present data should not be viewed as having broad applic- 
ability. They are restricted to a particular configuration 
of heater and container, and their value lies in that they 
illustrate the induced-convection effect. 

3. The expressions 


linax 

^maxp 

= f (W , I) 

(54a) 

^max 

1max F 

= f (W' , N) 

(55a) 


are the appropriate expressions to correlate q max data for 
the present heater configuration, in a larger container, 
over ranges of pressure, gravity, and size (as long as W'>1), 
and for different liquids. 

4. The peak heat flux for a horizontal ribbon heater (and prob- 
ably for other geometries as well) approaches a constant min- 
imum fraction of q ma x-p when the scale parameter, I and W' , 
become large. It also approaches this limit as the induced- 
convection buoyancy parameter, N, increases. 


66 


VII. SUMMARY OF DESIGN RESULTS 


Assessment of the RoTe of Gravity 

Let us now return to the initial motivation for this study — 
namely that of determining the interacting influences of gravity 
and configuration on the extreme heat fluxes. The preceding three 
chapters have shown beyond any doubt that 


%ax ^min 

or 

^max-p Iminp 


4>(L T , I) 


(57) 


is the appropriate expression 7 for correlating q max and q~ in . 

The use of equation (57) assumes that pressures are not close to 
the critical pressure; that L' is on the order of magnitude of 
unity or greater' and that surfaces are clean. 


Equation (57) is an adaptation of equation (20) . For the 
present discussion it will be chosen in preference to equation 
(16) which uses the induced buoyancy parameter, N. Thus the g- 
dependence of the expression is limited to qmaxp and L', and our 
subsequent manipulations will be significantly simplified. Re- 
writing equation (57) as 


^max or ^min “ ^max F or ^minp) $ 
we can form the following derivative: 

^ (imax or q m j_ n ) (qmaxp or Iminp) $ 

dg g 


1 1 d&n <j> 

-J- 

4 2 din L* 


(57a) 


(58) 


But 


^maxp or < 3minF )< f >=(< 2max or <3min> * so equation (58) becomes 


dJln (q max or qmj n ) _ 3. 1 d&n <j> 

din g 4 dJln L 1 


(59) 


Equation (59) displays the influence of g on the extreme 
heat fluxes in a very evocative way. It shows that the basic in- 
fluence of gravity is to increase q max or q m in as g^/^ (i.e. the 
lead term on the right hand side is 1/4) . However there is a 
second term — usually negative in sign--which can strongly detract 
from the gl/4 influence. This is the d£n$/2d5,nL' term which re- 


. 7 We introduce the special symbol <j> at this point to denote the 
ratio (qmax/qmaxp) or (gmin/^minp) / instead of the general func- 
tion, f. 
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fleets the additional influence of the scale parameter, L' . 


The transition between the ranges in which gravity augments 
the extreme heat flux, and detracts from it is also of major in- 
terest. This can be found by setting the right hand side of equa- 
tion (59) equal to zero, whence 


din 4> _ 1 

din L* ~ 2 


(60) 


Accordingly we have the following criterion: 

d ^n 4> r> “1/2, q ma x or q m i n increases with g T (60 * 
din L' 1 < -1/2, q max or q m in decreases with g J ' 

Thus, to satisfy our interest in the influence of gravity 
we should prefer to see the collected design results of the pres- 
ent work presented on in<J> versus inL' coordinates. We shall as- 
semble the results in this way after first exploring one more 
geometry that we have not yet mentioned. 


<3max and gmin on Spheres 


No sphere data have been obtained in the present study; 
however, a limited number of data obtained by other investigators 
provide the basis for further extending our set of correlation 
surfaces. Merte and Clark [6]; Lewis, Merte and Clark [7]; and 
Korayem [43] give a total of nine usable q m i n data for nitrogen 
at atmospheric pressure. Reference [7] also provides about 14 
qmax points for spheres . 

The q m -i n data were assembled in reference [3] under a 
slightly different scheme of correlation. We have correlated 
the same data on <j> versus R' coordinates in Fig. 31, and we have 
correlated the qmax data in the same way in Fig. 32. 

The resulting correlation lines must be accepted hesitantly. 
The data are restricted to nitrogen at one atmosphere . Thus , if 
any systematic deviations of nigrogen data occur within the _ pres- 
ent scheme of correlation, they will be entirely reflected in 
these lines. Furthermore there are too few data to clearly de- 
fine the behavior in the small R' regime. 

In the large R' limit the data rather strongly suggest that 


limit q ma x 
R’ 03 qmax F 


0.84, 


and 


limit q m jn 

R' ^ 00 qmin F 


0.45 


(61) 
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denotes uncertainty of point 


dimensionless radius, R’ 


Fig. 31 Minimum heat flux on spheres in nitrogen at I atm. 
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Fig. 32 Peak heat flux on spheres in nitrogen at I atm. 




for spheres. The limit of hydrodynamic behavior, or the R' be- 
low which the present scheme of correlation no longer applies, 
is not specified. However, we doubt that it should be valid very 
much below a diameter equal to half of a wavelength, or below 
2.73R'. 


Collected Dimensionless Heat Flux Curves 

Figure 33 presents the collected results for spheres and 
cylinders on tf> versus L 1 coordinates. Some selected values for 
q^ax on ribbons for I c =3200 are included for comparison on this 
plot. Let us note some of the important general characteristics 
of these curves : 

a. ) All curves are monotonically decreasing in L' . This 
means that, to the best of our knowledge, the "quarter-power" de- 
pendence of the extreme heat flux on g is an upper limit. Most 
configurations result in a slower increase of extreme heat flux. 

b . ) The great success of Zuber 1 s equation for qmaxp i n de- 
scribing q m ax f° r a variety of configurations is borne out in 
these curves. In all cases that we have considered q m ax is with- 
in about 50 percent of q m axp* 

c. ) As L'-*- 00 , the quarter-power dependence is re-established. 
This fact has not been borne out with data in the case of q m i n on- 
cylinders, since no such data exist for sufficiently large values 
of R' . 

d. ) q m i n exhibits much greater deviations from flat plate 
behavior than qmax does . 

e. ) qmin f° r cylinders actually decreases with g in the 
range that we have investigated. 
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I. = 3200 


0.07 0.1 0.2 0.4 12 4 

dimensionless radius, R 1 (or W 1 for ribbons) 


Fig. 33 Collected dimensionless extreme heat flux results. 





VIII. CONCLUSIONS 


The present studies were originally motivated by the need to 
understand how gravity and size interact in their influence on the 
peak and minimum pool boiling heat fluxes. A sequence of prior 
studies [10, 2, 3] had verified that q max and <3min could be corre- 
lated by 

gSgji- and qmin = f(L') (8) 

g maxp ^min-p 

under a wide variety of conditions. The dimensionless character- 
istic length, L', is 

L' = L[g(p g -p f )/CT] 1/2 (7a) 

This idea had served to explain the limited success of Zu- 
ber's flat plate equations in variable gravity as noted by Cos- 
tello, Siegel, Clark and others. Earlier investigators had at- 
tempted to compare their variable gravity data with the Vg depen- 
dence in Zuber’s equation, but without understanding the depen- 
dence of q max upon f (/g) that enters through the scale parameter, 
L', the effect of gravity could not be completely accounted for. 

Our first efforts under the support of NASA have therefore 
been directed specifically at a critical study of equation (8) . 

A large portion of these efforts are reported in references [11, 
12, 13, 14, and 15], and are reviewed in the preceding Chapters. 
They have briefly led to the following findings: 

A. Dimensional analysis combined with consideration of the 
physical nature of the problem leads to a general correlation 
which can be written in the form: 


?max ^min _ ^ , T , 

rf or J= : = f ( L , 

^irLinp 


I or N, /1+pg/Pf , 6 C ) 


(16 or 20) 


The roles of 0 C and of /l+p g /pf are conjectural in this expression. 
We believe that they can exert an influence under appropriate cir- 
cumstances, but we have not provided documentation of such influ- 
ences in the present studies . 


B. In the studies of film bo iling, there was no liquid- 
solid contact, the role of /1+p /pf was presumably incorporated 
in a nd the liquid flow pattern probably did not induce dis- 

ruptive 11 convective currents. Thus equation (16 or 20) reduced to 
equation (8) in this work. 


■ The wavelength of vapor removal during film boiling was 

found to vary as 


73 


<300 


o-l/2 

X d /ld F = (1+1/2 [ (l+b/R)R'] 2 ) 

where R' is an L' based upon the wire radius, A d p is the most sus- 
ceptible wavelength on a flat plate, and the vapor blanket thick- 
ness, b, is usually a minor factor which can be estimated by the 
method of Baumeister and Hamill. A variability of +60 percent and 
-25 percent in A was found to be inherent by virtue of the almost 
neutrally-stable character of the Taylor instability. Equation 
(30) was compared with data for many fluids over large ranges of 
g and R. 

The "best estimate" of qmin was found to be 

<3min = q minp t 0 . 0217/R* 2 (2R' 2 +l) ] (35) 

where q m i nF is Berenson's minimum heat flux for a flat plate, 
equation (14) . 

Zuber's "wave-collapse" mechanism of vapor removal was found 
to give way to a bubble-merger mechanism that was totally dominated 
by capillary forces, over the range: 0 . 07^R' £0 . 12 . When R' ex- 

ceeded unity, the "wave-collapse" mechanism deteriorated, although 
wavelengths equal to A dF could still be identified up to R'^2.5 — 
the upper limit of the observations . 


C. The peak heat flux on horizontal cylinders has been pre- 
dicted analytically. The result is in the form of equation (8) 
since I, 6 C , and (1+pf/pg) exert negligible influence: 


and 


Qfrnax _ 6 (R 1 +A) 2 /2 

%axp tt 2 /3 R ' 


0.15 <1 R ' £ 3.47 


^max 

^maxp 


0.894, R* > 3.47 


(43) 


where A is an L' based upon the vapor blanket thickness. The fol 
lowing empirical expression for the dimensionless blanket thick- 
ness was written to represent data over the range 0 . 15^R' £3 . 47 : 


A = [2.54R' + 6.48R' exp (-3 .AA/bT) ] 2 / 3 -R' 


(49) 


Combining equations (43) and (49) gave the following expression 
which converges almost to 0.894, and is valid for all R'>0.15: 


- maX = 0.89 + 2.27 exp (-3 . 44 /r 1 ’) (50) 

q max F 

These results conform with the simple version of our correlation 
strategy — equation (8) . 
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More than 400 original data points, and a comparable number 
cJf data from a variety of other sources verified equation (50) 
within ± 20 percent. 

D. Application of equation (16 or 20) in the correlation of 
about 900 original data for horizontal ribbons on an insulating 
base revealed a strong influence of X or N, especially for large 
N's or small I's. This influence was illustrated graphically in 
Figs. 27, 28, and 29. These Figures reinforced another finding 

of the present studies — namely that as the scale parameters, L' 
and I , become large , qmax/<3max F approach constant values . 

E. One of the chief overall conclusions to be drawn from 
the preceding studies is that an extraordinarily successful gen- 
eral scheme of correlation has been developed. However, Figure 
21 shows that it will not serve to predict the extreme heat 
fluxes when R' becomes small. For cylinders the method fails when 
R'<0.15. Our horizontal ribbon data only reached down to W 1 / 2-. 35, 
but at that point q max was also beginning to rise sharply and 
scatter broadly. Thus we are tempted to conclude that in any ge^ 
ometry the present correlation will fail when L'2 (the ratio of 
gravity to capillary forces, or "Bond number") falls below the 
order of 1/10 . 

Thus it remains for us to identify the correct limiting be- 
havior for q max under very low gravity. The correlations and pre- 
dictions that have been developed are useful in specifying q max 
and qxnin for systems that are not too small, under low gravity. 

But the basic change of behavior that occurs at small L' prevents 
extrapolation of results in that direction. 

F. As long as R' is not so small as to invalidate our 
method of correlation, the influence of gravity on qmax and qmin 
can be characterized by 

d&n (qmax or qmin) _ l i_ d 
dg “42 dL* 

The combined results of the present study, shown in Fig. 33, in- 
dicate that this derivative can become negative for only one of 
the cases that we have studied, namely qmin on horizontal cylin- 
ders . 


Jin 


qmax 

‘Jmaxt 


or 


qmin 


^mini 


(59) 
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APPENDIX A NOMENCLATURE 


4 H 


a 

b 

c v 
f (x) 

g 

G 


h fg 


t 


f g 


i 

k 

L 

L' 

N 

P r 


q 

^max 


vnaxp 


^min 


cross-sectional area of a vapor jet during strong nu- 
cleate boiling on a horizontal cylinder 

area of heater surface within the "unit cell" subtended 
by a vapor jet 

amplitude of an interfacial wave (see Fig. 6) 
thickness of the vapor blanket in film boiling 
specific heat at constant volume 
any arbitrary function of x 
acceleration of gravity 

acceleration of gravity under earth normal conditions 

g/g e 

latent heat of vaporization 
corrected hfg, equation (33) 

induced convection scale parameter, [pfLa] -*-/ 2 /y 

I based on W c 

wave number, 2t\/X 

any characteristic length 

L [g (pf-pg) /a] , dimensionless characteristic length 

induced convection buoyancy parameter, I 2 /L’ 

reduced pressure, pressure of system divided by critical 
pressure 

heat flux 

peak nucleate pool boiling heat flux 

qmax on an infinite flat plate (after pg. 1, q iaa x Tr is 
taken to be the value defined by equation (3) , whenever 
it is used) 

minimum film pool boiling heat flux 
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Sfitiirip 


R 

R' 



U 


gc 


W 

w c 

W' 

Y 


6 


A 

AT 


®c 

A 

X d 

x d F 

A 

A BH 

a lw 

’ a sk 

U 


q m in on an infinite flat plate (after pg. 12, qminp is 
taken to be the value defined by equation (14) whenever 
it is used) 

radius of a horizontal cylindrical heater 
Rig (pf-pg) /a] , dimensionless radius 

velocity of liquid approaching a heater 
velocity of vapor leaving a heater 

critical Ug for which liquid-vapor interface becomes 
Helmholtz unstable 

width horizontal ribbon heater 

width of test capsule 

W[g (pf-pg)/cr] x / 2 , dimensionless width 

dimensionless peak heat flux for R'<<1, (q max or q m i n )/ 
Pg h fg (cf/lj) 

effective vapor blanket thickness during strong nucleate 
boiling on a horizontal cylinder (see Fig. 14b) 

5 [9 (Pf”Pg) /cf] / dimensionless vapor blanket thickness 

heater surface temperature minus saturated liquid temp- 
erature 

contact angle among liquid, vapor, and the heater sv 
face 

wavelength 

most susceptible (or "most dangerous") Taylor unstable 
wavelength during film boiling 

A^ for a plane horizontal interface 

x d/ X d F 

Baumeister-Hamill A- based on equation (31) 

Lienhard-Wong A, based on equation (23) or (30) 
Siegel-Keshock A, based on equation (29) 
liquid viscosity 
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p f' p g 
a 

<t> 

0 ) 

a 

^max 


saturated liquid and vapor_densities , respectively 
surface tension between a liquid and its vapor 

^max/ ^maxp ) or ^min/ ^minp ' 
frequency of a wave in the complex plane 

dimensionless frequency defined by equation ( 26 ) 

Q for the maximum value of (iw) 
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